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Abstract 

A new version of the method of Continuum-Discretized Coupled-Channels 
method (CDCC) is proposed, that is, the Eikonal-CDCC method (E-CDCC). 
The E-CDCC equation, for Coulomb dissociation in particular, can easily and 
safely be solved, since it is a first-order differential equation and contains no huge 
angular momentum in contrast to the CDCC one. The scattering amplitude 
calculated by E-CDCC has a similar form to that by CDCC. Then one can 
construct hybrid amplitude in an intuitive way, i.e., CDCC amplitude is adopted 
for a smaller angular momentum L and E-CDCC one for a larger L related to an 
impact parameter b. The hybrid calculation is found to perfectly reproduce the 
quantum mechanical result for 58 Ni( 8 B, 7 Be+p) 58 Ni at 240 MeV, which shows its 
applicability to systematic analysis of 8 B dissociation to extract the astrophysical 
factor Sn with high accuracy. 

1 Introduction 

The solar neutrino problem is one of the central issues in the neutrino physics [1 . 
Nowadays, the neutrino oscillation is assumed to be the solution of the problem and 
the solar neutrino physics moves its focus on determining oscillation parameters: the 
mass difference among z/ e , and u T , and their mixing angles [2]. The astrophysical 
factor 5*17, defined by Si 7 (E) = cr pi (E)E exp[27rr/] with a pl the cross section of the 
p-capture reaction 7 Be(p, 7) 8 B and n the Sommerfeld parameter, plays an essential 
role in the parameter-search procedure, since the prediction value for the flux of 8 B 
neutrino, which is intensively being detected on the earth, is proportional to Si7(0) . 
The required accuracy from astrophysics is about 5% in errors. 

Because of difficulties of direct measurements for the p-capture reaction at very low 
energies, alternative indirect measurements were proposed, p-transfer reaction [HI EJ E] 
with the Asymptotic Normalization Coefficient (ANC) method |H] and 8 B Coulomb 
dissociation IE1 E] are typical examples of them; we concentrate on the latter in the 
present paper. 8 B Coulomb dissociation can be assumed to be the inverse reaction of 
7 Be(p, 7) 8 B provided the 8 B is dissociated through only its El transition by absorption 
of virtual photons. Once this condition is satisfied, cross section of the p-capture reac- 
tion, thus Sn(E), can easily be obtained by applying the principle of detailed balance 



to the measured dissociation cross sections. Intensive measurements of the 8 B Coulomb 
dissociation are being made in RIKEN jj], GSI [Sj and MSU 0; the extracted 5i 7 (0) 
are almost consistent each others. There is, however, contradiction between MSU and 
RIKEN data about the contribution of the E2 component. Additionally, role of nuclear 
interaction, interference between nuclear and Coulomb interactions in particular, has 
not yet clarified quantitatively. Therefore, it can be said that the extracted 617(0) from 
those 8 B Coulomb dissociation measurements contain uncertainties to some extent. 

In order to determine 5^(0) accurately, careful analysis of the 8 B dissociation in- 
cluding both nuclear and Coulomb interactions is necessary. The method of Continuum- 
Discretized Coupled- Channels (CDCC) [TQJ , which was proposed and developed by 
Kyushu group, is suitable for that purpose. In fact, it was shown in Ref. f9j that 
CDCC can very well reproduce the MSU data. It is not straightforward, however, to 
extract the component corresponding to the inverse reaction of 7 Be(p, 7) 8 B out of the 
total measured spectra. 

Very recently [IT] , we proposed a procedure that determines 5i7(0) from the 8 B dis- 
sociation measurements by using the ANC method, which is free from uncertainties by 
the use of the detailed balance. An important advantage of the CDCC + ANC analysis 
is that one can quantitatively evaluate the accuracy of the extracted 5*17(0); the fluctua- 
tion of ANC by changing the 8 B single particle wave functions is interpreted as the error 
of 5*17(0) coming from the use of the ANC method. We analyzed 58 Ni( 8 B, 7 Be+p) 58 Ni 
at 25.8 MeV measured at Notre Dame [Ej, for which detailed balance was found to fail 
to obtain 5i7(0) because of its low incident energy [THj. The extracted 5i7(0) is 22.83 
± 0.51 (theo) ± 2.28 (expt) eVb and the ANC method turned out to work very well, 
i.e., less than 1% of error; the remaining theoretical error comes from the choice of the 
modelspace of CDCC calculation. Although quite large systematic error (10%) of the 
experimental data prevents one from determining 5^(0) with the required accuracy, 
this method significantly reduces theoretical errors of 5^(0) . 

The CDCC + ANC analysis is expected to be applicable to the RIKEN and MSU 
data at several tens of MeV/nucleon. From practical point of view, however, CDCC 
calculation including long-ranged Coulomb coupling-potentials requires extremely large 
modelspace rather difficult to handle; typically the required number of partial waves is 
15,000 for the MSU data jHJ- Although interpolation technique for angular momentum 
reduces the number of Coupled-Channels (CC) equations to be solved, CC equation 
with huge angular momentum L is rather unstable and careful treatment is necessary. 
In this sense, it seems almost impossible to apply CDCC to the GSI data at 250 
MeV/nucleon, where more than 100,000 partial waves will be required. 

In the present paper we propose a new version of CDCC; the Eikonal-CDCC method 
(E-CDCC) being based on eikonal and non-adiabatic calculation. Since E-CDCC equa- 
tions contain no angular momentum for the relative motion of a projectile and a target 
nucleus, the 5 matrix is obtained easily and safely. Another important feature of the 
method is that the resultant scattering amplitude has the form which is very similar to 
that obtained by CDCC, by using the relation between impact parameter b and L. This 
allows one hybrid calculation, i.e., one can adopt the amplitude of E-CDCC only for 
larger L, where classical picture is expected to be valid, and connect them with that of 
the standard, fully quantum mechanical, CDCC for smaller L. The hybrid calculation 



includes all quantum-mechanical effects necessary through the CDCC amplitudes. In 
addition to that, the use of the E-CDCC amplitudes removes all problems concerning 
with huge angular momenta and drastically reduces computation time. From theo- 
retical point of view, the hybrid calculation will give an insight into the connection 
between quantum-mechanical and classical pictures. 

In Sec. ["""formalism of E-CDCC and the construction of hybrid scattering amplitude 
are briefly described. We show in Sec. [""] the validity of the hybrid calculation for 
58 Ni( 8 B, 7 Be+p) 58 Ni at 240 MeV; results for the elastic and total breakup cross sections 
are compared with those of fully quantum-mechanical calculation. Effects of adiabatic 
approximation and importance of the interference between the larger- and smaller- L 
regions are also discussed. Finally summary and conclusions are given in Sec. EJ 

2 Formalism 

In this section we briefly describe formalism of E-CDCC and how to construct hybrid 
scattering amplitude by using a result of (standard) CDCC partly. Detailed formalism 
and theoretical foundation of CDCC are shown elsewhere [""T"| 111*1 Ho]. 



Figure 1: Schematic illustration of the system concerned in the present paper. 

The system that we treat in the present study is illustrated in Fig. 1. We start with 
the expansion of the total wave function 



where t is the total spin of c and m is its projection on the z-axis taken to be parallel 
to the incident beam; the subscript represents the initial state. $j^ m is the channel 
wave function of c, which contains both bound and scattering states and for the latter 
discretization of the continuum is made. {Qi t e m } is assumed to form an approximate 
complete set for a finite configuration space being significant for a reaction concerned. 
We here make the following eikonal approximation: 
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where c denotes channels {i, £, m} together and the wave number K c is defined by 



with b the impact parameter; the direction of K c is assumed to be parallel to the z-axis. 

Inserting Eqs. (JTJ) and (J2J) into the three-body Schrodinger equation neglecting the 
second order derivative of ifj c , one can obtain 

5-X?>-%?>(*) = i>f(Ae <{Kf (4) 

where J-f£) is the form factor, i.e., the interaction between A and each constituent of c, 
folded by <3> c and $ c '; we put b in a superscript since it is not a dynamical variable but an 
input parameter. Equation (J3J is solved with the boundary condition oo) = <5 c0 . 

Since the E-CDCC equation is first-order differential one and contains no coefficient 
with huge angular momentum, one can easily and safely solve the equation. 

Using the solution of the E-CDCC equation, the scattering amplitude with E-CDCC 
is given by 

where /i is the reduced mass of the c + A system. Making use of the following forward- 
scattering approximation: 

e *(KW-KiW).R = - K ^sme } bcoscj> R +{Kf -K^cosd f )z 

« -K®6 f bcos<t> R + {Kf -K$*)z, (6) 

one obtains 

(7) 

where the eikonal S'-matrix elements are defined by Sf^ miQiQmo = ^^(oo). 
We here discretize f E : 
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where bf n , 6£ ax and 6£ id are defined through xf^bf 11 = L, K^hf™ = L + 1 and 
Kiim = L + l/2, respectively. In deriving Eq. (|SJ we neglecting the 6-dependence 
of KtfJn, exp[— iK^ m 6fb cos 0^] and <S^j ^ omo within small size bin of b corresponding 
to each L. After manipulation one can obtain 



f E ~ 2n V K ° J 2L + 1 i (m - mo) Y, (K')(s {bfid) -6-6**6 ) 



(9) 



which has a similar form to that of standard CDCC: 



/3mM>mo = JfT Yl £ S \j L °^ 1 (L o 0£ o m o \Jm o )(L m -m £m\Jm ) 
tK o l j=\L-e\ L a =\j-£ \ v 47r 

x (Sf UM - ^ALoWR m ~ mo ^„o(K'). (10) 
The construction of the hybrid scattering amplitude / H is rather straightforward: 

fi£m,i £ m = Yl II + T"! ff.l (H) 
L<L C L>Lc 

where (ff) is the L-component of (/ E ). -^c represents the connecting point be- 
tween quantum-mechanical and classical pictures, which is chosen so that coincides 
with for L > Lq. It should be noted that Eq. (|TT|) includes all quantum-mechanical 
effects being necessary and also interference between the two L-regions. 

In the above formulation we neglect the Coulomb distortion. In order to include it, 
we use xi{R-> @r) ~ Xc^Pi z )0c(^j z ) instead of Eq. (J2J, where 0^ is the Coulomb wave 
function The formulation of f E can then be done just in the same way above. 



3 Results and discussion 

In order to see the validity of the hybrid calculation with CDCC and E-CDCC, we 
analyze 58 Ni( 8 B, 7 Be+p) 58 Ni at 240 MeV. Parameters of the modelspace are as follows. 
The numbers of bin-states of 8 B are 16, 8 and 8 for s-, p- and d-states, respectively. The 
maximum excitation energy of 8 B is 10 MeV and r max (-R max ) is 100 fm (500 fm). L max is 
taken to be 1000 for saving computation time, which is somewhat small compared with 
the expected value of L max that gives perfect convergence, i.e., L max ~ 4000. However, 
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Figure 2: Angular distribution of the elastic (left panel) and total breakup (right panel) 
cross sections for 58 Ni( 8 B, 7 Be+p) 58 Ni at 240 MeV. The solid, dashed and dotted lines 
show, respectively, the results of quantum-mechanical (QM), eikonal (EK) and hybrid 
(HY) calculation with Lq = 100. The dash-dotted line in the right panel represents 
HY calcualtion with Lq = 200. 



the validity of the hybrid calculation can definitely be clarified in the modelspace 



above. As for 8 B wave functions, the single-particle model by Kim et al. was 
adopted. For nuclear interaction between 7 Be (p) and 58 Ni we used the parameter set 
of Cook et al. [TH] (Becchetti and Greenlees jUj); we neglected the spin-dependent part 
for simplicity. Additionally, we replaced Eq. (jSJ) by (h 2 /2fi)K 2 = E — e iy i and neglected 
c-dependence of K in the calculation of f E , which was found to have no effects on 
numerical results. We stress that this is not an adiabatic approximation since K c is 
explicitly treated in solving the E-CDCC equation (HJ). 

In the left and right panels in Fig. 2 we show the elastic cross section (Rutherford 
ratio) and total breakup cross section as a function of scattering angle in center-of mass 
(cm.) frame, respectively. The solid, dashed and dotted lines represent the results of 
quantum-mechanical (QM), eikonal (EK) and hybrid (HY) calculation, where the Lq 
are taken to be 1000, and 100. In the right panel the result of HY calculation with 
Lq = 200 is also shown by the dash-dotted line. The agreement between QM and 
HY calculation with appropriate value of Lq, namely 100 (200) for elastic (breakup) 
cross section, is excellent; the error is only less than 1%. One also sees fairy large 
difference between EK and QM calculation. Since our EK calculation contains no 
corrections to the straight-line approximation, this does not directly show the fail of 
EK approximation. However, it seems quite difficult for EK calculation to obtain 
"perfect" agreement with the result of QM one. On the contrary, the HY calculation 
turned out to be applicable to analyses of 8 B dissociation where very high accuracy is 
required. 

HY(Lc = 200) 
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Figure 3: The p-states breakup 
cross sections by QM calculation with 
(dashed line) and without (solid line) 
adiabatic (AD) approximation. The 
maximum number of L is taken to be 
2500. For comparison non- adiabatic 
QM result with L max = 1000 is shown 
by the dotted line. 



Figure 4: The p-state breakup cross 
sections by QM calculation with < 
L < 200 (dashed line) and by EK cal- 
culation with 200 < L < 1000 (dotted 
line). The dash-dotted line is the in- 
coherent sum of the dashed and dotted 
lines and the solid line is the HY result 
with Lq = 200, namely, the coherent 
sum of the two. 



We show in Fig. 3 the p-state breakup cross sections by QM calculation with (dashed 



line) and without (solid line) adiabatic (AD) approximation. Since the convergence of 
AD calculation is rather slow compared with non-AD one, we take L max = 2500 to make 
fair comparison. One sees that AD calculation much overestimates the breakup cross 
section, which shows the failure of AD approximation for this reaction. In Fig. 3 the 
non-AD QM result with L max = 1000 is also shown by the dotted line. The difference 
between the solid and dotted lines is finite but quite small, which indicates the QM 
calculation with L max = 2500 works well, i.e., the large difference between AD and 
non-AD results is indeed due to the use of AD approximation. 

In Fig. 4 the p-state breakup cross sections by QM calculation with < L < 200 
and EK calculation with 200 < L < 1000 are shown by the dashed and dotted lines, 
respectively. The dash-dotted line is the incoherent sum of the two, which deviates 
from the HY result shown by the solid line. This result shows that the essence of our 
HY calculation is the construction of the hybrid scattering amplitude not the hybrid 
cross section. 

4 Summary and Conclusions 

In the present paper we propose a new method to treat breakup process accurately 
and efficiently, that is, the hybrid (HY) calculation with the Continuum-Discretized 
Coupled-Channels method (CDCC) and the Eikonal-CDCC method (E-CDCC). E- 
CDCC describes the center-of-mass (cm.) motion between the projectile and the 
target nucleus by straight-line (eikonal approximation) and treats the excitation of 
the projectile explicitly, by constructing discretized-continuum states as in CDCC, i.e., 
non-adiabatic and non-perturbative calculation can be done. When Coulomb distortion 
appears, we make a similar approximation to the eikonal one; we use Coulomb wave 
functions instead of plane waves. The resultant scattering amplitude by E-CDCC has 
a similar form to the quantum-mechanical (QM) one obtained by CDCC, which allows 
one to construct HY amplitude / H in an intuitive way; / H is given by the sum of the 
partial amplitude calculated by CDCC for smaller angular momentum L and that by 
E-CDCC for larger L, i.e., larger impact parameter b. 

We analyzed 58 Ni( 8 B, 7 Be+p) 58 Ni at 240 MeV by QM, full-eikonal (EK) and HY 
calculation. It was found that the elastic and total breakup cross sections obtained by 
QM calculation are "perfectly" reproduced by the HY calculation, namely, the error 
is only less than 1%. Constructing HY cross section, not HY amplitude, turned out 
to fail to reproduce the corresponding QM result, which shows the importance of the 
interference between the two L-regions above. 

In conclusion, the hybrid calculation using CDCC and E-CDCC allows one very 
accurate analyses for breakup processes. The accuracy of the model is enough to be 
applied to 8 B dissociation relating to the astrophysical factor 617(0), the aim of which is 
determining S^O) with less than 5% errors. E-CDCC drastically reduces computaion 
time and eliminates many problems concerned with huge angular momentum in solving 
Coupled-Channels (CC) equations. Thus, the hybrid calculation opened the door to the 
systematic analyses of 8 B dissociation measured at RIKEN, MSU and GSI. Extracted 
617(0) by using the ANC method will be reported in near future. 
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